Phase diagram for a single flexible magnetic filament in a poor solvent 
at zero field 
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The equilibrium conformations of a flexible permanent magnetic chain that consists of a sequence of linked magnetic colloidal 
-nanoparticles with short-ranged attractive interactions are thoroughly analysed via Langevin dynamics simulations. A tentative 
phase diagram is presented in the limit of infinite dilution. We find that the phase diagram of such colloidal polymers exhibits 
several unusual conformational phases when compared with the non-magnetic chains. These are characterised by a large degree 
]of conformational anisotropy: simple closed chains, helicoidal-like states, partially collapsed states, and very compact disordered 
states. The phase diagram shows several interesting features like the existence of at least two 'triple points'. 



1 Introduction 

Artificial magnetic filaments can be obtained by mutually 
linking magnetic colloids to form a chain. These magnetic 
chains represent the equivalent to magnetic polymers but at 
supra-molecular scale. Magnetic polymers have magnetic 
properties only at T<100Kii^, magnetic filaments can retain 
their magnetism at zero field and room temperature if the size 
of the nanocolloids is chosen adequately. Although the ex- 
.act range may depend strongly on the type of materials used, 
common magnetic cores with sizes between 10 and lOOnm 
can exhibit a permanent dipolar moment if they are single 
"magnetic domains blocked in the crystal lattice of the parti- 
cle^. Cores smaller than lOnm are expected to exhibit a super- 
paramagnetic behaviour due to the dominance of Neel relax- 
.ation processes. Cores larger than lOO^mare expected to show 
paramagnetic properties due to the formation of multiple mag- 
netic domains^. 

The path towards the synthesis of such permanent magnetic 
filaments has been possible thanks to a progressive increase in 
the abilities to control the size of magnetic colloids and the na- 
.ture of bonds between the colloidal particles ^ . Initially, Furst 
"et al.^ reported in their seminal papers the use of micron- 
sized magnetic-filled paramagnetic latex beads as monomers 
to form magnetic chains. Further progress was achieved by 
Dreyfus et al. ^ using DNA to link micron- sized paramagnetic 
beads. The first steps towards the assembly of magnetic fila- 
jnents using sub-micrometric particles (0.5 — 0.8/im) where 
reported by Tabata et al.^ Magnetic chains with lengths in 
the range 30 to 50/im, either in solution or attached to a sur- 
face, were presented by Singh et al. ^^'^^ and Toonder et al. 
Shorter magnetic chains synthesised using sub-micrometric 
colloids in an electrolyte solution ^^-^^ and via single- step syn- 
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thesisi^ have also been achieved. On the other hand, the syn- 
thesis of arrays made of long magnetic filaments (from 1 to 
200/im) with length-to-monomer- size ratios up to 125 times 
were obtained by Evans and co-workers 

Further technological improvements allowed the use of 
nanoparticles with sizes below 50^m that form mesoscopic 
one-dimensional structures—"—. Very recently, Benkoski et 

24,25 j^^^g reported the optical visualisation of small fer- 
romagnetic cobalt nanoparticles (23.5^m) organised into mi- 
croscopic one-dimensional chains without needing to resort 
to templates or multiple fabrication steps. Those filaments 
can be attached to a surface and form brush-like structures by 
conveniently arranging some magnets along the surface2i2^. 
Very recently, Sarkar and Mandal^^ have performed a note- 
worthy synthesis of magnetic chains using DNA as a template 
on which they directly grow the magnetic nanoparticles with 
sizes ranging between 7 and 17 nanometers. The particles re- 
main anchored to the DNA, and the whole system behaves as a 
magnetic filament. It is also worthwhile to point out that there 
have been successful attempts by Zhou et al. to lock and pre- 
serve the structural conformations of filaments made of mag- 
netic cobalt nanocolloids of 20^m in size^^. Goubault^^ et 
al. have achieved the synthesis of flexible magnetic filaments 
by the simple procedure of bridging the surfactant layers that 
ferrofluid particles carry adsorbed on top of their surface in 
presence of an external magnetic field. Once the bridging has 
occurred, the filaments are irreversibly formed and the exter- 
nal magnetic field can be removed. 

The emerging interest in this relatively novel field is due 
to the fact that magnetic filaments are very appealing from 
the technological point of view. They can be thought as im- 
proved substitutes of current ferrofluids, or as elements for 
magnetic memories, chemical and pressure nanosensors, or 
medical applications, to mention just a few^. Furthermore, 
the study of the structures adopted by magnetic filaments is 
also of relevance in the area of inverse ferrofluids where non- 
magnetic particles are added to a ferrofluid solution. Those 
particles will behave inside the ferrofluid equivalently to that 
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an ordinary magnetic particle would do in a non-magnetic sol- 
vent By linking a sequence of non-magnetic nanocolloids 
to form chains and immersing them in a ferrofluid, we can 
have exactly the same behaviour for those colloidal chains as 
for true magnetic filaments in an ordinary non-magnetic sol- 
vent. The use of such 'polymer colloids' in ferrofluids can 
have also many potential applications. For instance, 'inverse 
filaments' have already been used in the assembly of photonic 
crystals that only exist during the time that an external mag- 
netic field is applied to the sample In general, in most of 
the applications, the knowledge of the different types of equi- 
librium structures that filaments may adopt can be of extreme 
importance. Nonetheless, so far, very little is known about 
the structures that permanent magnetic chains may adopt as a 
function of temperature, length, and other related parameters 
like magnetic momentum of the particles, or the strength of 
the attractive interactions among colloidal particles. 

In the case of non-magnetic attractive chains the study of 
their phase transitions has been exhaustive, see refs. ^i-^'^, and 
refs. ^^'^^ for a review. In the case of semiflexible attractive 
polymers it is known that there exist several conformational 
phases that are different from the typical swollen coils, col- 
lapsed globules, crystal and glassy states. Several woks have 
reported semiflexible chains to exhibit toroidal or disk-like 
phases Helix structures have also been found for some 

very specific square- well potentials'^. Related to it, helical 
long-lived transient states have also been identified for chains 
with truncated Lennard- Jones potentials'^. 

For magnetic chains the number of studies is much lower. 
In addition to the previously mentioned works devoted to 
the experimental synthesis of magnetic filaments, there have 
been several attempts to obtain phase diagrams for magnetic 
chains using Ising or Heisenberg-like monomers in good sol- 
vent—"—, corresponding to the case of non- attractive colloids. 
Henceforth, we will refer to such kind of chains as non- sticky 
filaments. In addition, some works have also dealt theoret- 
ically with the investigation of the magneto statics of chains 
made of magnetic nanoparticles of different shapes The 
derivation of the partition function, the intra-chain correla- 
tions, and the coil-globule transition for flexible non- sticky 
magnetic chains in the limits of zero and infinitely strong ex- 
ternal magnetic fields has also been pursued A large frac- 
tion of the existing works has been devoted to the study of 
the properties of magnetic filaments made of paramagnetic or 
super-paramagnetic non-sticky chains. Most of these stud- 
ies treat the filaments as elastic rods to use them as micro- 
propellers (micros wimmers)'-^^'^'"^^ under the action of an 
external field, or as actuators ^^'^' with the purpose of perform- 
ing tasks similar to those of micrometric magnetic cilia^^*^, 
but at the nanoscale. 

For lower dimensionalities, Sanchez et al. have recently 
addressed via numerical simulations the adsorption and the 



equilibrium conformations that result when a non- sticky mag- 
netic filament is adsorbed on top of an attractive surface. In 
their work a preliminary phase diagram was sketched for those 
filaments as a function of the temperature and the strength of 
the surface attractive interaction. In the adsorbed regime two 
different states were observed: at low temperatures, closed- 
chain structures were found. Such closed-chains are more or 
less extended chains with the two ends at close contact. At 
larger temperatures, the chains showed open conformations 
that resemble those observed for ordinary non-magnetic ad- 
sorbed chains. An even more recent work from the same au- 
thors has characterised the different structural regimes dis- 
played by a non- sticky flexible magnetic filament immersed 
in a perfect solvent when varying the strength of the dipolar 
coupling parameter. 

Other works related to the study of the structures adopted 
by single magnetic filaments, specially at low temperatures, 
are those devoted to the study of clusters of free particles 
that interact via Stockmayer potentials (Lennard- Jones plus 
point dipole potentials). Thus, for instance. Miller et al.^^ 
have found that clusters of particles interacting via Stock- 
mayer potentials exhibit a rich variety of ground states that 
includes rings and different types of coils with several topo- 
logical knots. The type of ground state exhibit by a cluster 
of Stockmayer particles was found to depend strongly on the 
dipole moment and the number of particles. Nonetheless, to 
our best knowledge, it is unknown if the conformations of sin- 
gle magnetic sticky filaments at low temperatures will shear 
some resemblance with the ground states found for Stock- 
mayer clusters which lack of any permanent links between 
particles. 

In this paper we investigate the influence of magnetic in- 
teractions on the phase diagram of sticky and non- sticky 
magnetic chains in three-dimensions via the use of extensive 
Langevin dynamics simulations. In Section [21 we describe the 
numerical model and the details of the simulations. In Section 
Owe present and discuss our results which, in sake of clarity, 
we have split in three main parts. The first and second part are 
devoted to the discussion of the structures and phases that exist 
in the weak and the strong regime for the attractive interaction, 
respectively. In the third part, all the results are collected in or- 
der to build up a tentative phase diagram. Finally, a summary 
and a discussion of the conclusions are presented in Section IH 

2 Numerical Model 

In order to investigate the behaviour of flexible magnetic fila- 
ments in the limit of infinite dilution we have performed nu- 
merical simulations in which a single magnetic filament is 
placed in a unbounded open space. The magnetic filament 
is modelled as a bead- spring chain made of a sequence of 
magnetic beads (colloidal particles) of diameter cr^, carrying 
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a point dipole Sit their centre. Henceforth, the subindex e 
denotes the experimental values of the physical quantities we 
use. The absence of such subindex means the quantity is ex- 
pressed in reduced units. We will use through all this work 
reduced units. Thus, for a length 4 its corresponding reduced 
value is I = le/Oe. In the simulations, in sake of simplicity, the 
diameter of the colloidal particles is set to cr = 1 , so all length 
scales are measured in units of this diameter. 

Two different types of particles will be considered: non 
sticky particles, where the predominant interaction between 
particles is the steric repulsion due to the their cores; and 
sticky particles, in which in addition to the repulsion there 
exists an attractive pair- wise interaction among the particles. 
Henceforth we will refer to this last interaction as 'the LJ in- 
teraction'. The attractive interaction between two particles / 
and j will be modelled via the following potential that com- 
bines the core repulsive part (cutoff rcut = and the at- 
tractive part (rcut = 2.5cj), 

Uatt (r) = VtsLj{r, (T, 1 , rcut = 2 ^ cj) +y,,L/(r, (J, e, rcut = 2.5cj) 

(1) 

where r is the distance between the centres of the particles 
/ and 7, i.e. r = \ri — rj\, and Vtsu is a truncated- shifted 
Lennard- Jones potential—. 
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where Uu{r) = 4£[(cr/r)^^ — (cr/r)^]. The LJ energy param- 
eter e, is given in units of the experimental well depth e^, and 
any energy U will be also referred to e^, i.e. U = Ue/Se- In 
the same way we choose the Boltzmann constant to be ^ = 1 
in reduced units, and therefore the reduced temperature is 
T = keTe/ £e- The modulus of the dipole moments can be also 
expressed in the reduced system as: ji^ = ji^ / {Anjio^e^^^^e)' It 
should be noted that the soft-core and the attractive part have 
been implemented as in eq. ([T]) and not through a simple LJ 
potential because we want to ensure that the effective repul- 
sion is roughly the same when different values of £ < 1 are 
used. In this way a comparison between chains with differ- 
ent depths for the attractive well and non- sticky chains can be 
performed more easily. 

The colloidal particles are assumed to interact pair-wise as 
point dipoles according to the potential. 
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where Tij = Vi — rj is the displacement vector between parti- 
cles / and j. The energy due to magnetic interactions is calcu- 
lated by direct summation over all pairs of particles. In spite 
of being algorithmically ^{N^), for small numbers of colloids 
this is the fastest and most accurate way to compute it due to 



the open boundary conditions. The values of /i = to be 
used depend in general on the composition and size of the col- 
loidal particles: aside from cobalt nanoparticles, colloidal par- 
ticles found in common commercial ferrofluids usually do not 
exceed values of /i 10. 

In order to connect the colloids to form a chain, a linking 
model for the particles slightly more involved than the one 
used in previous works ^^'^^ has been implemented. In this 
new approach springs between consecutive particles are not 
anchored at the centres of the beads, but at fixed points on 
their surface, as shown schematically in Figure [T] The rea- 
son to implement a novel variant of the bead- spring model 
is to better mimic the real behaviour of magnetic filaments 
made of sticky colloidal particles. In our preliminary tests us- 
ing springs anchored at the centre of the dipoles, we noticed 
that, for moderate and strong magnetic coupling, the sticky 
filaments form packed structures holding filament sections in 
which the dipole orientation is inverted with respect to the 
chain backbone in order to minimize the overall magnetic en- 
ergy. Nonetheless, in most common real magnetic filaments, 
those inversions would imply a strong stretching of the bond 
linking the colloidal particles at the borders of the two differ- 
ent domains. Thus, the existence of such inverted domains 
in real filaments must be heavily penalised. Instead, chains 
modelled with springs anchored at fixed points on the surface 
of the beads, penalise such 180 degrees rotations in an effec- 
tive way. The model we present is therefore useful to analyse 
the behaviour of particles with permanent magnetic moment, 
i.e. those that having a magnetic monodomain that is blocked 
in the crystal structure of the particle and in consequence the 
dipole is only allowed to relax via the Brownian mechanism. 

The proposed spring potential is written as: 



(4) 



where and rj are the position vectors of the centres of the 
beads. Ui and Uj are unitary vectors placed along the direction 
that joins the two anchoring surface points of each bead (see 
figure [TJ. Thus, the anchoring points are collinear and located 
at I'l = UiG/2 and li = —UiO/l with respect to the centre of 
the bead. We assume all links in the chain are formed accord- 
ing to the following scheme: the point on the surface of the 
/ — 1 particle with position r^-i + is linked to the point on 
the surface of the particle / with position + U. In order to 
obtain the desired effect of severely penalising conformations 
with consecutive dipoles in anti-parallel configuration, we as- 
sociate each director vector Ui with the dipole moment of the 
particle, i.e. Ui = The constant of the potential 

is set to = 30 which is enough to ensure the average bond 
length to lie within a reasonable range r^^^j G [0.98, l.ljcr. 
The use of larger values for Ks is possible but implies a further 
reduction of the integration time step. 
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(a) (b) 



Fig. 1 The magnetic filament is modelled as a chain of beads linked 
by springs anchored onto a point of their surfaces, see eq. ©. The 
magnetic moment of each particle is used as reference vector to 
define the unitary vectors u = fi/jj. which are used to determine the 
position of the anchoring points on the surface of the particles, 
g/2 u and — (7/2 m, when the centre of the particle is taken as the 
origin. 



The numerical simulations are performed using Langevin 
dynamics, in which colloidal particles are moved according to 
the translational and rotational Langevin equations of motion 
that for a given particle / are^: 

Mi^=Fi-rTVi^^J (5) 

Ir^ = r,-rua)t^^f (6) 

where F^, and Ti are respectively the total force and torque 
acting on the particle /, Mf and If are its mass and iner- 
tia tensor while Ft and Tr are the translational and rota- 
tional friction constants. and are Gaussian random 
forces and torques, each of zero mean and satisfying the 
usual fluctuation-dissipation relations. In the simulations, 
t = te \J I {j^e<^}\ where rUe is the real mass of the colloids; 
F = FeGe/ and T = Te/^e- For equilibrium simulations, the 
values of the mass, the inertia tensor, as well as friction con- 
stants r^, and Tr are irrelevant because the same equilibrium 
state is reached independently of their value. Only the dy- 
namics to attain such equilibrium state may show differences. 
Thus, in sake of simplicity, the particle mass is chosen to be 
m = 1 and we take the inertia tensor to be the identity matrix 
in order to ensure isotropic rotations 7 = 1. We have chosen 
= 1 and Tr = ?>/A because we observed that these values 
produced a conveniently fast relaxation to equilibrium 
The reduced time step is set to 5r = 5 • 10~^ in order to en- 
sure a correct integration of the equations of motion. 

The simulations are started by placing the filament in a open 
three-dimensional non-bounded space with the position of the 
first bead located randomly. The remaining monomers are po- 
sitioned using a self-avoiding random walk scheme with an 
overlap radius of 0.9cr. The chain is pre-equilibrated at T = 1 
for 2-10^ integrations with the magnetic interaction turned 



off while the time step is slowly increased from 10~^5r till 
0.055r. Subsequently, magnetic interactions are turned on, 
and a second pre-equilibration stage consisting of 5 • 10^ in- 
tegrations is performed while gradually raising the time step 
from 0.1 • \0~^5t till 5^. Right after, if the final temperature 
is r < 1 we perform an annealing process using the final time 
step St: the temperature is reduced from T = 1 down to its fi- 
nal value by performing a set of five annealing stages of 5 • 10^ 
steps each one. Once the final temperature has been reached, 
the chain is equilibrated for a period of 3 • 10^ ^^/^ 5t m order 
to ensure that the chain is in the thermodynamic equilibrium. 
After the equilibration period, the system is sampled at inter- 
vals of 2500 e^/^ 8t for another period of 15 • 10^ e^l^ 5t to 
make sure there are no correlations between measurements. 
To get further assurance that the results do not depend on the 
initial conditions and to improve statistics, an average over 15 
independent runs for each set of sampled parameters (r, e, /i) 
is performed. The simulations have been performed using the 
package ESPResSaii. 

3 Results and discussion 

In the present model there are two main competing interac- 
tions: on the one hand there is the LJ attractive interactions 
among the beads that tend to collapse the chain when the tem- 
perature is lowered^2i2S; on the other hand there is the mag- 
netic interaction which is known from ferrofluid studies ^^'^^ to 
favour the formation of rod-chains and rings in which dipoles 
tend to align in a nose-tail conformation. Nose-tail conforma- 
tions are those in which two dipole / and j satisfy jU^ • flj = /i^, 
and Hi • rij = Mj ' ^ij = ^ij- This noise-tail alignment 
allows to minimise the magnetic energy in eq. ©. Thus, 
it seems reasonable to define the dimensionless parameter 

3 

T] = ^ which measures the relative strength of the attrac- 
tive LJ interaction compared to the strength of the magnetic 
interaction for particles at close contact, rij = O", in a nose-tail 
conformation. 

Although we just focus on the behaviour of a single chain, 
the number of parameters involved to explore the full phase 
diagram is large. Therefore, in order to get a first sketch on 
a section of the the phase diagram, we have chosen to focus 
on a specific case with a fixed number of beads N = 100. In 
addition, all beads are assumed to have the same magnetic mo- 
ment jl^ = 5 and share the same LJ energy parameter e. In this 
first approach no external magnetic field is present. The phase 
diagram will be studied as a function of two parameters: i], 
that is controlled modifying the value of e while keeping ji 
fixed, and the temperature T. The range of values explored 
for the parameter t] is T] G [0,0.2]. T] = corresponds to the 
case in which particles do only exhibit sterical repulsion be- 
tween their cores and no effective attractive LJ interaction ex- 
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ist among them (non- sticky filaments). This case is quite com- 
mon when magnetic particles are sterically stabilised against 
flocculation, as is the case of ferrofluids^^. The range of tem- 
peratures sampled is T G [0.27,5], where the upper boundary 
r = 5 has been chosen because the most interesting features 
were found to occur well below such value in all cases. The 
relative strength of the magnetic to the thermal energy, is usu- 
ally expressed in terms of the dipolar coupling parameter A, 
which is defined as half of the ratio of the magnetic energy 
of two dipoles in close contact and in a nose-tail alignment to 
the thermal energy, i.e. A = /i^/ {Trl^^^). Thus, the value of 
the effective dipolar coupling parameter for the systems under 
study ranges in A G (0.6, 20), being A 10 for T = 0.5, and 
A ^ 1 for r = 3.5. Nonetheless, even at T = 5 the thermal mo- 
tion is still not powerful enough as to completely overcome the 
magnetic interactions: the main structural observables of the 
chain have values still far from those corresponding to an ideal 
self avoiding random walk that is the behaviour one could ex- 
pect in the limit T ^ oo. Nonetheless, once the chain is in 
a extended open state, we do not expect any other relevant 
change to happen by raising T. 

Values of T] > 0.2 and T < 0.27 correspond to regions in the 
phase space characterised by compact structures. An effective 
sampling of these regions with usual Langevin methods is very 
costly in computer terms and require special techniques like 
e.g. umbrella sampling Wang-Landau method ^^'^^ or other 
advanced existing methods and is left for future work, 
although our preliminary results suggest the existence of in- 
teresting features that are worthwhile to be explored 

In sake of a better understanding we have split the discus- 
sion of the results in three parts. In the first part we will discuss 
the behaviour of the magnetic filaments for T] < 0.10, in the 
second part for T] > 0. 10 and in the third and last part we will 
discuss the implications our results have for the phase diagram 
of magnetic filaments. 

3.1 Filaments with t] < 0.10 

The radius of gyration Rg and the end-to-end distance Ree are 
two important observables that can be very useful in order to 
follow the structural changes of a magnetic filament. The end- 
to-end distance is defined as Ree = {{n — r^q)^)^/^ where (...) 
denotes an average over all the sampled conformations of the 
chain. On the other hand, we can define the gyration tensor 
through their elements, 

where a and j3 denotes de Cartesian components x, y, and 
z. The tensor can be represented as a diagonalisable 3x3 
matrix with three eigenvalues or principal moments hence- 



forth labelled as Af > A^ > A3 . The radius of gyration is 

7?^ = ^Af + A| + A|. 

Figure [2] shows the end-to-end distance (top plot) and the 
gyration radius (bottom plot) as a function of the temperature. 
Our results show that the ends of the chain tend to get closer 
as the temperature is lowered. Remarkably, the behaviour of 
the end-to-end distance is very similar for all filaments with 
7] G [0,0.07]. In the range 7] > 0.07 - 0.10 a noticeable two- 
fold decay step emerges whose origin is explained below. 

Also, as can be observed in figure O the radius of gyra- 
tion of sticky filaments (r] > 0) shows a progressive differ- 
ence, when compared with non-sticky filaments (t] = 0). For 
all measured values of t] > 0, the magnetic filaments tend to 
shrink when there is an attractive LJ interaction among the 
beads. As in the case of the end-to-end distance, the radius of 
gyration also tends to develop a two-fold decay for t] 0.10. 

In the case of non-sticky particles, i] = 0, the variation in 
the values of Rg is quite small. It ranges from a value of 
Rg ^ 9g 8itT = 5 (value not shown in the figure [2]) to a value 
of Rg ^ 12(7 at r ~ 0.3. Nonetheless the behaviour of Rg 
is clearly different from the one expected for non-magnetic 
chains: for such chains without attractive interactions, one 
would expect Rg to be approximately constant if the bond 
length is fixed. Here, as the temperature is lowered, the mag- 
netic filaments show an initial expansion, followed by a con- 
traction T e [0.7, 1.0] and a second expansion at T < 0.7. This 
particular behaviour is also observed in filaments with values 
of T] closer to T] = down to a temperature in which the at- 
tractive interactions dominate and induce a strong collapse of 
the chain. This expansion-contraction behaviour can be un- 
derstood as follows: as the temperature is lowered the mag- 
netic interactions for an extended open chain dominate, and 
favour the stretching of the chain in a conformation in which 
all the dipoles tend to remain aligned along the same direction. 
Nonetheless, below a certain temperature, the most favourable 
conformation is a closed structure because an extra aligned 
pair can be created by getting closer the ends of the chain. The 
entropic penalty of a closed chain is not excessive when the 
temperature is low enough. The closed structures at low tem- 
peratures are already known to occur in ferrofluids where 
particles assemble into ring-like structures. Although the ten- 
dency is to form a ring when the temperature is lowered, the 
structure of the filament is far from being a perfect ring at the 
point where the filament closes. This can be observed in the 
plots corresponding to T] = in figure [3] An additional proof 
in favour of this relation between the closing of the chain and 
the slight decay of the radius of gyration is the fact that both 
phenomena occur in the same range of temperatures (see fig- 
ure [2]). For low values of T] the transition from open to closed 
structures occurs at T ~ 0.85, which corresponds to a value of 
the dipolar coupling parameter of A ~ 6. These results are in 
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Fig. 2 The plot depicts the the end-to-end distance Ree (top) and the 
radius of gyration Rg (bottom) as a function of the temperature for 
?7 < 0.10. The case r\ =0.10 has been also included for a better 
comparison with the results corresponding to T] > 0.10. 



good agreement with what has been observed for non-sticky 
chains in the simpler variant of the bead- spring model 

Further insight on the behaviour of the magnetic filaments 
in the regime r\ < 0. 10 can be obtained by examining the spe- 
cific heat Cy, defined in terms of our dimensionless tempera- 
ture as 



Cv 



(8) 



where U is the total energy of the filament [U = Uatt + U^np + 
Us\ . The experimental specific heat is related to this definition 
according to Cye =keCy. Figure m shows the specific heat as 
a function of temperature. The bottom plot in figured shows 
some small peaks appearing around T 0.8 — 0.9 which can 
be identified with the transition from open to closed structures. 
Three main features can be observed for such peaks: the first 
one is that for values of T] G [0,0.07] they almost coincide in 
position, value and width. This means that the corresponding 
transitions are almost independent of the strength of the at- 
tractive interaction. This independence is due to the absence 
of close contacts between particles which are not first near- 
est neighbours. The second feature is that compared to the 
peaks that can be observed in the top plot of figure lU the peaks 
around T ~ 0.8 — 0.9 represent very small fluctuations in the 
energy. This fact is coherent with the idea that those peaks 
represent a transition from extended open to simple closed 
structures because the difference in the total energy between 
one state and the other should not be much larger than the en- 
ergy due to the creation of a new pair of aligned dipoles, plus 
the energy due to the close contact of two particles. A third 



feature is that when the value of r] increases, the peaks repre- 
senting the transition from open to closed structures become 
partially overlapped by the tails of the larger peaks present in 
the specific heat. This overlap makes it very difficult to de- 
termine the temperature of the open-closed transition from the 
peaks of the specific heat for values of t] > 0.07. In fact, the 
only observable we have found to roughly determine the point 
of this transition is dRee{T)/dT. It is reasonable to consider 
that the transition approximately coincides with the inflection 
point of the end-to-end distance, i.e. dRee{T) / dT has an ex- 
trema at that temperature. Nonetheless the scenario becomes 
more complex due to the fact that, as rf ^0.10, Ree shows a 
two-fold decay. We will come back to this point later in sec- 
tion 13.31 where a global picture of the phase diagram will be 
available. 

The large peaks observed in the specific heat in figure |4] 
represent a different type of transition, namely, the conversion 
of the simple closed structures into compact helicoidal-like 
structures as it can be observed in figure |3] for snapshots at 
rf =0.05 and T < 0.6. As it is shown in figure [5] the main part 
of the filament adopts a structure that remembers a tight he- 
lix, while the ends of the chain arrange in such a way that the 
two ends stay in close contact. These helicoidal states are re- 
lated to the toroidal conformations observed in non-magnetic 
semi-flexible chains where local chain stiffness helps to sta- 
bilise those structures ^^'^^ . In our case, the magnetic interac- 
tions tend to force a nose-tail orientation of the dipoles that, 
in addition, will induce a local chain stiffness Nonetheless, 
there are some subtle differences between the toroidal con- 
formations found in non-magnetic chains and the helicoidal 
structures observed here: in the non-magnetic case the toroidal 
walls are thick with a width of several particle diameters while 
in the magnetic case they tend to be much thinner. An open 
question is what would happen if longer chains of the order 
of A^- 10^- 10"^ were studied. We argue that in the case of 
magnetic chains, an helicoidal state is preferred to a toroidal 
conformation, at least if 7] is not very large. An helicoidal 
structure allows to minimise the energy associated to pairs of 
dipoles with their dipole moments lying parallel = = M 
but for which ri2 • II = 0. Unlike the nose- tail conforma- 
tions, these pairs have the highest possible magnetic energy 
and therefore are heavily penalised. Such energy penalty de- 
creases as ~ with the distance between the two dipoles. 
In a toroidal conformation, a similar number of unfavourable 
dipole pairs may exist, but in difference to the helicoidal struc- 
ture the distance between the two particles will be, in general, 
much shorter, leading to a higher energy penalty. The situation 
reverts for large values of T] in which the short ranged attrac- 
tive interaction dominates over the magnetic one and a torus 
is preferred. Thus, helicoidal structures seem to be the result 
of a complex interplay between the attractive interactions that 
tend to collapse the chain into an isotropic globule, the mag- 
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7^=0.00, T=0.300 7^=0.00, T=0.475 7^ = 0.00, T=0.600 7^ = 0.00, T=0.900 7^ = 0.00, T=1.100 




7^=0.05, T=0.300 7^=0.05, T=0.475 7^ = 0.05, T=0.600 7^ = 0.05, T=0.900 7^ = 0.05, T=1.100 




Fig. 3 The figure depicts several typical snapshots of the configurations that magnetic filaments of length N =100 adopt for different values 
of the temperature and relative strength of the attractive interactions rj. Snapshots for T] = 0, 0.05 are shown. The two ends of the filaments 
are show in green and red colours. The rest of the beads are painted in two colours (yellow and blue) to show the orientation of the magnetic 
moments of the particles. 



netic forces which, on the one hand, tend to locally stretch 
the chain by leading to an effective local stiffness that favours 
toroidal conformations but, on the other hand, tend to avoid 
the formation of pairs of parallel dipoles with their relative 
vector position perpendicular to the direction of the dipoles, 
and the chain entropy. It is also worth to mention that heli- 
coidal states have been found for non-magnetic chains when 
specific short-ranged square- well potentials were used^^. An- 
other issue worth to mention is the possibility that those heli- 
coidal states could be long-lived metastable transient states as 
those found by Sabeur et al. for simple homopolymers using 
truncated Lennard- Jones potentials whose monomers where 
subject to simple spring potentials anchored on their centres. 
Since we use similar reduced units, and share the same simu- 
lation technique (Langevin dynamics), a rough comparison is 
feasible: Sabeur et al. observed that the decay from the helical 
states is a stochastic rate-driven process, where the escape rate 
is l/to ^ exp{AE/kT) and AE is the height of the energy bar- 
rier. In their simulations they found a value of to ^ 2000 5^ for 
the particular case of a homopolymer chain of length N = 100 
at r = 0.04. Since our equilibration and measurement times 
are of the order of 3 • 10^ e^^^ 5t and 15 • 10^ e^^^ 5t respec- 
tively, and taking into account that the lowest temperature we 
have sampled is of the order of T ~ 0.25, we can reasonably 
conclude that the helicoidal structures we observe are not ex- 
pected to be long-lived metastable states but true equilibrium 



states. 

Figured] shows that the peaks in the specific heat signalling 
the formation of the helicoidal states are much larger than 
those related to the extended open to simple closed structure 
transitions. This fact is a consequence of the substantial in- 
crease in the number of particles at close contact that appear in 
the helicoidal state with respect to the simple closed structures, 
in which the main interactions are only between the nearest 
neighbours in the chain sequence and between the chain ends. 
The global chain rearrangements show up as large fluctuations 
in the total energy, and thus, larger peaks in the specific heat 
are expected. Another feature of the helicoidal transition is 
that the transition temperatures in this case show a strong de- 
pendence with T] , which is quite opposite to the behaviour ob- 
served for the transition from extended open to simple closed 
structures where the transition temperature remains almost in- 
variant. 

It is worth to remark that, within the range of sampled 
temperatures, peaks associated to such transition are only 
observed for T] > 0.02. We expect such helicoidal transi- 
tions to take place also at t] < 0.02 and smaller temperatures 
(r < 0.27). A question that remains open is whether the tran- 
sition from simple closed to helicoidal conformations takes 
place for all non-zero values of t] , or if there is a critical value 
for the strength of the attractive interaction, T]^, such that for 
7] < r]c the closed-helicoidal transition is not possible. 
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A quantitative way to characterise the formation of an he- 
Hcoidal state is to evaluate a typical order parameter able to 
signal the formation of an helix. Among the different order pa- 
rameters^^, we have chosen the so-called H^^ parameter which 
characterises the global helical order defined as 

^4 = T7^( lV.- - ^.--i) X (r,-+i - n)). (9) 

^ i=2 

^4 = is associated to isotropic conformations that resemble 
a rod, whereas H4 = I holds for perfect helix. Other order 
parameters like were found to lead to similar conclusions 
that those derived from H4. 

In figure [6] (top plot) the value of the order parameter H4 as 
a function of the temperature is depicted for t] < 0.10. The 
values of H4 remain equal to zero until temperatures close to 
the helicoidal transition point. The derivative, dH4/dT, shows 
that the position of the inflection points occurring at the high- 
est temperature coincides with the position of the peaks in the 
specific heat. This fact shows that both observables are linked 
to the appearance of helicoidal states. Figure |6] shows that 
the achievement of the helicoidal states is very gradual: the 
further the temperature is lowered the higher is the value of 
H4, and the structure looks more similar to a perfect helix. 
Only for values of T] ^ 0.10, the order parameter H4 seems 
to reach a plateau within the range of temperatures sampled. 
In any case, the largest values of H4 ~ 0.4 are quite low com- 
pared with those of an ideal helix H4 = I. These relatively 
small values for the order parameter have a two-fold cause: 
the first one is that only the main part of the chain adopts an 
helix-like structure, whereas the ending parts of the chain ar- 
range in such a way that the ends can be at close contact. If 
the whole structure was in an helicoidal state the two ends 
would be separated by a large distance, which gives an energy 
penalty. Such small differences in energy may be irrelevant 
at high temperatures but not at low temperatures. The second 
reason lies in the fact that these helicoidal states do not look 
like normal cylinders but rather they exhibit a symmetry sim- 
ilar to that of an elliptic cylinder, whereas the parameter H4 
is defined having in mind ideal helixes with a symmetry sim- 
ilar to that of a regular cylinder. Notice that the cross product 
{rt — Ti-i) X (r/+i — Ti) would be zero for a very elongated 
ellipse in which bonds between particles are locally aligned. 

The asymmetry of the filaments can be roughly inferred 
from a visual inspection of the different conformations shown 
in figure [3] However, a good set of observables to determine 
quantitatively the degree of asymmetry of the different confor- 
mations are the ratios of the second and third eigenvalues of 
the gyration tensor to the main eigenvalue, / and A3 / Af . 
Those ratios are shown in|7]as a function of temperature. In the 
region of high temperatures T e [1.8, 5] the ratios of the eigen- 
values are approximately constant and equal to A^/A^^ ~ 0.2 
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Fig. 4 The specific heat as defined in eq. lU is represented as a 
function of temperature for different values of the relative strength 
of the LJ attractive interactions T] . The bottom plot is a zoom to 
highlight the small peaks that exist at T ~ 0.8 — 0.9 which point out 
the transition from extended open chains to simple closed structures. 




(a) (b) 



Fig. 5 The figure depicts two different perspectives of a typical 
conformation in the helicoidal state (t] = 0.10, and T = 0.475). The 
code of colours follows from figure |3] 
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Fig. 6 The helicoidal order parameter H4 as defined in eq. ^ is 
plotted as a function of temperature for several values of T] < 0.10 
(top plot) and T] > 0.10 (bottom plot). 



Fig. 7 The ratio of the second and third eigenvalues of the gyration 
tensor, see eq. 0, to the main eigenvalue are shown as a function of 
the temperature for several values of r] < 0.10. 



and A^/A^^ ~ 0.06, respectively. These ratios mean that fil- 
aments are highly elongated along one direction and almost 
lie in a plane because the value of third eigenvalue is almost 
negligible respect the first and second eigenvalues. These be- 
haviour roughly corresponds to the tendency of the magnetic 
particles to align themselves in a row. The fact that the value 
of the ratios is almost the same for all r] < 0. 10 implies that in 
the range of high temperatures the attractive interaction plays 
a very minor role. Nonetheless, we are still far from the limit 
T ^ 00 vvere we expect to observe filaments with values of 
the eigenvectors that resemble those of an ideal SAW with 

Another interesting feature is that the eigenvalue ratios for 
< 7] < 0.07 are almost identical to the case r] = up to 
the helicoidal transition takes place. This implies that in the 
transition from simple closed to helicoidal states, the shape 
of the closed structures for sticky filaments with T] G (0,0.07) 
is similar to that of simple closed structures formed by non- 
sticky filaments at the same temperature. In the region where 
helicoidal structures exist, a comparison between the sticky 
filaments with the non- sticky ones (t] = 0) shows that the at- 
tractive interaction, as it could be expected, favours the more 
isotropic structures. 

A comparison of figures |7] and |2] shows that the temperature 
at which the eigenvalue ratios separate from the non-sticky 
case coincides with the temperature at which Rg suffers a more 
prominent decay. Thus, when the chain collapses by reduc- 
ing its size in every direction, the global shape becomes more 
symmetric. 

In the region where compact helicoidal structures exist, the 



second eigenvalue exhibits a substantial increase, A2 / A^ ~ 
0.74, which shows that the helicoidal states do not project on 
top of a circle but rather of an ellipse. This fact is consistent 
for all the values of 17 sampled. Below a certain temperature, 
the degree of anisotropy of the shapes of the filaments grad- 
ually increases again, a signal that points in the direction that 
another conformational phase is arising. We shall come back 
to this point in section [T2l 

At first sight, it can be surprising that all closed and col- 
lapsed structures do not show, in average, an isotropic distribu- 
tion of particles, at least along their two main axis of rotation. 
This can be understood as follows: in addition to an ideal dis- 
tribution of N particles over a circumference with their dipole 
moments tangent to the curve, there exist other closed struc- 
tures which are non isotropic and have an energy that is com- 
parable to that of an ideal circumference. One should expect 
the system at finite T to sample different states of similar en- 
ergy, and consequently, the averaged values of the eigenval- 
ues of the gyration tensor should also reflect a certain degree 
of anisotropy. In addition we should take into account that 
bonds between particles are not constant, but fluctuate around 
their mean value which can further enhance the number of 
anisotropic states with energies similar to those of a ring con- 
formation. 

3.2 Filaments with 77 > 0. 10 

As we have already shown in figure lU the specific heat devel- 
ops an additional peak for t] = 0.1 at T ~ 0.4. This peak is 
also observed for t] > 0. 1 (see top of figure[8]) and, as r] ^ 0.2, 
it shifts towards higher temperatures T G [0.3,0.7] while de- 
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Fig. 8 The figure depicts in the top plot the specific heat as a 
function of the temperature for several values of the relative strength 
of the attractive interaction t] > 0.10. The bottom plot shows the 
derivative of the gyration radius dRg/dT as a function of the 
temperature for the same range of T7's . 




Fig. 9 The radius of gyration (top plot) and the end-to-end distance 
(bottom plot) are depicted as a function of temperature for several 
values of the relative strength of the attractive to the magnetic 
interaction r\. 



creases in height. The radius of gyration and the end-to-end 
distance depicted in figure [9] show that those peaks in the Cy 
must correspond to a transition that takes place when the chain 
is already in a collapsed state. For that reason it is very dif- 
ficult to discern from the structural parameters any relevant 
signal of the transition. This transition must involve an inter- 
nal rearrangement of the particles without noticeable changes 
in the global size of the structure. A rough idea of such inter- 
nal rearrangements is provided by the helical order parameter 
^4 depicted in the bottom plot of figure [6] for 7] > 0.10. A 
comparison between Cy and //4 (top of figure [S]) reveals that 
the inflection point of the lines in the interval T G [0.3, 0.7] 
roughly coincides with the position on the peak in the specific 
heat. We can infer that the new transition is associated with 
a loss of the helicoidal order and the onset of compact disor- 
dered states. A typical structure in this regime is represented 
in figure [TOk. 

The knowledge of the eigenvalues allows to make a rough 
estimation of the volume that the magnetic filament occupies 
if an assumption about the shape of the filament is made. For 
the most compact case studied, that corresponds to 7] =0.2 
in the low temperature regime {T G [0.3,0.6]), we found the 
eigenvalues of the gyration tensor to be approximately 4, 
Al ~ 2, and = 1. The ratios of the eigenvalues as well as 
the snapshots of the systems at such temperatures suggest that 
plausible shape for the filament is an elliptic cylinder. Under 
this assumption, the major semi-axis is a ^ vT^) + 1, the mi- 
nor semi-axis is /? ^ + 1 and the height is /z ^ VT^) + 1 . 



The correcting factor + 1 is added because the gyration tensor 
only takes into account the centre of the particles, and the cal- 
culated elliptic cylinder will neglect about half of the volume 
of the most outer particles. Once the parameters of the elliptic 
cylinder are known, the calculation of the approximate volume 
of the elliptic cylinder is straightforward: = nabh ^ 1 14cr^. 
To get an idea of the degree of compaction achieved, we 
should compare the previous volume to that of a compact 
hexagonal packaging (hep) ofN = 100 spheres, Vhcp ^ 74(7^. 
This result shows that the magnetic filaments at 7] =0.2 and 
r < 0.6 exhibit a degree of compaction close to that of a hep 
lattice of impenetrable spheres. On the other hand, for r] > 0. 1 
we observe that the helicoidal structures developed are far 
more isotropic than for T] < 0.1. Thus, for T] = 0.2 the highest 
ratios of the eigenvalues are A^/A^^ ~ 0.9 and A3 /A^^ ~ 0.6 
that imply a significant increase in the level of isotropy when 
compared to those values for t] < 0. 1 that are plotted in figure 

El 

From the comparison among Cy, Rg and Ree for t] > 0.1 
we can extract more useful information about the magnetic 
chains. A remarkable feature in figure |9]is that Rg and Ree 
have a very similar dependence with the temperature. This 
result differs from the observed behaviour for t] < 0.10 in 
which the filament first suffers a closing transition as T is low- 
ered and only then, by further reducing the temperature, there 
is a collapse of the chain towards a compact helicoidal state, 
/. e. the end-to-end distance decays at higher temperatures than 
Rg does. However, for 17 > 0.10, the attractive interaction is 
strong enough to force compaction to occur at higher temper- 
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atures than those at which filaments would suffer the closing 
transition. In principle, the closing temperature should not be 
substantially modified respect to the value for the non sticky 
case because the value of the magnetic dipoles is kept con- 
stant. Nonetheless, once the compaction of the chain occurs, 
the distance between the chain ends must substantially de- 
crease, and this lets the closing transition to occur right after 
the compaction. 

Another remarkable feature that can be extracted from fig- 
ure |9] is the observation that for the curves corresponding 
to r/ = 0.15 and rj = 0.20 it is possible to clearly asso- 
ciate the inflection points in Rg{T) with the two emerging 
peaks in the specific heat (figure [8] top) at T ~ 1.3 — 1.4 and 
T ^ 1.8 — 1.9, respectively. This behaviour suggests the exis- 
tence of a regime of intermediate states between the expanded 
chains and the compact helicoidal configurations. For values 
of 7] G [0.10,0.15) one can infer that the peak is also present 
but hidden in the long tail associated to the appearance of he- 
licoidal states that take place in the interval T G [0.8, 1.1]. For 
these smaller values of T] the fingerprint of the transition to- 
wards this intermediate regime is found in the double stage de- 
cay of both Ree and Rg (figure[9]). A natural question that arises 
is what kind of conformations do exist in such intermediate re- 
gion. The plots of Rg and Ree in figure [9] show that for those 
intermediate states, the chain is still far from being fully col- 
lapsed and there is still, on average, a long distance between 
the two ends of the chain. Snapshots of such intermediate 
states, as the one shown in figure [TOb, confirms the previous 
suggestion: between expanded open chains and the compact 
helicoidal states there exists a region of partially collapsed fil- 
aments in which a part of the filament is already in a compact 
state while the other parts are still in a expanded conformation, 
with the chain ends being separated by a relatively long dis- 
tance. Those partially collapsed structures remember vaguely 
of the core-shell structures found in non-magnetic semiflex- 
ible attractive chains Nonetheless in the present case, we 
have a core plus some loose tails rather than a shell surround- 
ing a core. This behaviour comes from the anisotropic nature 
of the magnetic interactions which favour relative straight seg- 
ments far from the core rather than the wrapping of the core 
by the non-collapsed part of the chain. It is not yet clear what 
the behaviourof the filament will be in the limit N ^ oo when 
long non-collapsed segments may exist. 

As we have shown above, the characterisation of this new 
transition from open expanded chains to partially collapsed 
states by using the specific heat is hard to be accomplished 
for T] < 0.20. We found that the best observable to deter- 
mine the transition point are the dRg/dT curves (see bottom 
of figure [8]). The maxima of the peaks of such function for 
T] = 0.15, 0.20, coincide with the apparent position of the 
emerging peaks for the Cy at the highest temperature. It is 
worth to mention that this observable has been also success- 




(a) (b) 

Fig. 10 The figure depicts two typical snapshots for high strengths 
of the LJ attractive interaction when compared to the strength of the 
magnetic interaction: (a) rf = 0.15, T = 0.40. (b) rf = 0.20, 
r = 1.60. The code of colours follows from figure |3] 



fully used in the determination of conformational phase transi- 
tions of non-magnetic chain s . Sharma and Kumar have 
verified that the locations of peaks in structural and thermody- 
namic quantities accompanying conformational transitions of 
macromolecules do coincide in the limit of A/^ ^ oo. For finite 
sizes, the two different types of observables may yield slightly 
different transition temperatures. However, in our case, the 
possible mismatches that could exist with respect the true tran- 
sition temperature are below our accuracy in determining the 
transition points, and we do not expect such slight mismatches 
to substantially modify our sketch of the phase diagram for 
magnetic filaments. Results for the transition temperatures 
obtained via dRg/dT are shown in figure [TT] as solid blue 
triangles-up for all values of rj . 

3.3 The phase diagram for isolated flexible magnetic fil- 
aments 

Gathering together the results presented in sections 13.11 and 
13.21 it is possible to build up a tentative sketch of the (T,ri)- 
phase diagram as shown in figure [TT] The solid black circles 
correspond to the transition points for t] < 0.10 derived from 
the position of the maxima of the peaks in the Cy and cor- 
respond to the transition from extended open chains to sim- 
ple closed structures (see bottom plot in figure |4]l. Solid red 
squares depict the transition points obtained from the maxi- 
mum of the highest peaks of the Cy in the range T G [0.7, 1 .2] , 
see figures 13] (top) and [8] (top). Those large peaks correspond 
to transitions towards a compact helicoidal state when the tem- 
perature is lowered. Solid green diamonds depict the tran- 
sition points from compact helicoidal states to compact disor- 
dered states, which are obtained from the maxima of the peaks 
of the Cy in the region of very low temperatures T < 0.7 in 
figure [5] (top). The solid blue triangles-up correspond to the 
maxima of the peaks in the dRg/dT, for T] > 0.1, that also 
mark the inflection points of Rg, (figure [8] bottom). As we 
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described in section [T2l the position of such peaks should ba- 
sically coincide with the transition points from open extended 
chains to partially collapsed states. Figure [TT] also shows that 
dRg/dT does a very good estimation of the transition points 
from simple closed chains to helicoidal states, in which tran- 
sition temperatures derived from dRg/dT are very similar to 
those obtained from the position of the peaks in the Cy . 

In figure [TT] the solid black circles only refer to the transi- 
tion for T] < 0.7. This is due to the fact that, for higher values 
of T] , the peak in the specific heat associated to such transition 
is hidden by the tail of the larger peak associated to the transi- 
tion to helicoidal structures. In section [TTI we mentioned the 
possibility of using the dRee/dT in order to characterise the 
transition from extended open chains to simple closed states. 
The inflection points of the Ree{T) curves are plotted in figure 
[TTIas solid magenta down triangles. Furthermore, for values 
of 7] < 0. 1 the inflection point in the end-to-end distance is 
clearly related to the transition point from open structures to 
simple closed structures. For values of T] > 0.1 the inflection 
in the Ree takes place approximately at the same temperature 
as the inflection point of Rg (solid blue triangles up) which, 
as discussed in section 13.11 is a consequence of the fact that 
a partial compaction of the chain triggers the closing of the 
chain. 

A very remarkable fact observed in figure [TTlis the existence 
of two different conformational 'triple points'. In the first 
'triple point' extended open chains would coexist with par- 
tially collapsed states and simple closed states. In the second 
'triple point', simple closed states will coexist with compact 
helicoidal states and partially collapsed states. The existence 
of two different, yet close triple points is a vivid example of 
how rich and complex the phase diagram for magnetic chains 
can be. 

Finally, it is interesting to discuss how the phase diagram 
would change if the strength of the dipolar interaction /i^ is 
modified, or if the range of the attractive forces is changed. To 
this end it is important to note that the bottom left corner of the 
phase diagram is dominated by the magnetic interactions, the 
upper left corner of the phase diagram is dominated by the LJ- 
like attractive interactions, and the bottom right comer is dom- 
inated by thermal motion. By increasing the dipolar strength 
11^ while keeping the other factors unaltered we should ex- 
pect the region where the magnetic interactions are dominant 
to expand. That means that the transition from extended open 
chains to simple closed structures should shift towards higher 
temperatures, and the transition from simple closed structures 
to compact helicoidal states should occur at higher values of 
rj . The region where partially collapsed structures exist should 
also shift towards regions of higher values of rj and T. A re- 
verse behaviour should be observed if we decrease the value 
of /i^ rather than increasing it. On the other hand, if we in- 
crease the range of the attractive interactions, the region where 
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Fig. 11 A tentative phase diagram for magnetic filaments of length 
N = 100 and /i^ = 5 is presented. See text in section [331 for a 
detailed explanation of how the different transition lines were 
obtained. In sake of clarity of the origin of each transition line, they 
have been painted using different colours and symbols. 



the attractive interactions dominate should expand, and there- 
fore one should expect the transition from compact disordered 
structures to helicoidal states to occur at lower values of r] . In 
turn, the transition from helicoidal to simple closed structures 
should also happen at lower values of 7]. It should be also 
possible to observe partially collapsed states at lower values 
of 7]. 

4 Conclusions 

A tentative phase diagram in the limit of infinite dilution 
for magnetic filaments made of N = 100 colloidal particles 
of identical size and magnetic momentum /i is presented as 
a function of the temperature and the relative strength be- 
tween the attractive LJ-like interactions versus the dipolar 
magnetic moment, T] . We have explored the range of parame- 
ters T G (0.27, 5) and T] G [0, 0.2] via numerical simulations 
using Langevin dynamics. Our results evidence a rich phase 
diagram in which it has been possible to characterise up to 
five different conformational phases and two 'triple points'. In 
the first triple point, located in the nearby region to (T ~ 0.9, 
T] ~ 0.1), open chains, partially collapsed states and simple 
closed states will coexist, while in the second 'triple point', lo- 
cated in the nearby region to (T ~ 0.6 -0.8,7]-- 0.07 - 0.09), 
the three phases in coexistence are simple closed chains, com- 
pact helicoidal states, and partially collapsed states. 

The characterisation of the different transition lines has 
been performed mainly via the analysis of the behaviour of 
the specific heat Cy, the radius of gyration Rg, and the end-to- 
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end distance Ree • At low temperatures three different phases 
are found depending on the relative intensity of the attrac- 
tive to the magnetic interaction: for low values of 7] sim- 
ple chain-closed structures are observed, whereas for inter- 
mediate values of rj compact helicoidal structures are devel- 
oped. Both phases exhibit a considerable degree of anisotropy 
along the two major axis, where the structure is much more 
elongated along one of them. The helicoidal states become 
more isotropic as T] increases. In the limit of strong attrac- 
tive interactions (large T] values), the helicoidal structure is 
lost and very compact disordered states with volumes close 
to the hexagonal-compact-package limit become the preferred 
states. Nonetheless, as figure [TOh illustrates, the filaments do 
not develop, in the range of sampled temperatures, a crystal 
order. 

Moreover at higher temperatures (T >l) and low 7] the pre- 
dominant states are extended open chains which exhibit a sig- 
nificant anisotropy due to the effects arising from the mag- 
netic interactions that favour the formation of flat structures 
with long persistence lengths. In these structures, dipoles are 
aligned in a nose-tail orientation. Nonetheless, as temperature 
increases, in the limit T ^ oo, the structures are expected to 
become more isotropic and resemble a swollen self-avoiding 
random coil. By increasing 7], a new kind of partially col- 
lapsed structure emerges. These structures are partially col- 
lapsed while keeping the two chain ends in a swollen confor- 
mation. The transition from extended open chains to partially 
collapsed states can be roughly characterised via the inflection 
points in the radius of gyration and the end-to-end distance. If 
the strength of the attractive interaction is further increased, a 
second transition from partially collapsed states to helicoidal 
states occurs. In this case, the transition can be monitored via 
the specific heat. This transition involves a severe rearrange- 
ment of the monomers and the peak associated to such transi- 
tion can be easily identified. Though not observed in the range 
of 77's sampled, we speculate that at even higher strengths of 
the attractive interaction, a new transition from helicoidal state 
to compact disordered globules should occur. 

Although the present phase diagram is only a simple sketch 
of a much more complex reality, interesting open questions 
emerge from it. One of them is, for instance, whether or not it 
is possible to find a critical 7] below which in the limit T ^0 
the only ground states are simple closed structures as those 
found for non- sticky filaments. The behaviour of the heli- 
coidal states for 17 0.2 is another issue worth to be stud- 
ied. As pointed out in section 13.11 helixes are energetically 
favoured over toroidal conformations for moderate values of 
7] , but for increasing values of this parameter there is the possi- 
bility that helixes convert into toroidal structures quite similar 
to those observed in semi-flexible non-magnetic chains ^^'^^ . It 
is yet unclear if compact helicoidal states will transit to com- 
pact globules directly, or rather they will do it via interme- 



diate partially collapsed states where the existence of a third 
'triple point' cannot be discarded. The structure of the ground 
states for the magnetic filaments, and up to which degree they 
show resemblances with the ground states observed for clus- 
ters of Stockmayer-particles^^, is also another open question. 
Finally, it will be also interesting to study how the phase di- 
agram changes with the chain length N or with the polydis- 
persity in size of the colloidal particles that form the magnetic 
filaments. 

A good understanding of the different structures that mag- 
netic filaments may adopt as a function of the interplay of 
the different interactions involved, and the phase diagram they 
display, is crucial in order to assess the use of these filaments 
for new technological applications or as substitutes of current 
ferrofluids with enhanced properties. Magnetic filaments have 
an enormous potential for new applications, and the charac- 
terisation of their properties is an open challenge. The present 
work represents a first step towards the understanding of the 
magnetic filaments that should stimulate further developments 
on this subject of increasing scientific interest. 
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